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ABSTRACT 


We  consider  the  inventory  control  of  a  single  item  in  a  two  echelon 
system  in  which  a  warehouse  (the  upper  echelon)  supports  N(N  _>  1)  re¬ 
tailers  (the  lower  echelon).  Customers  return  units  in  a  repairable 
state  as  well  as  demand  units  in  a  serviceable  state  at  the  retailer  level 
only.  The  return  and  demand  processes  are  assumed  to  be  mutually  indepen¬ 
dent.  The  stationary  return  rate  is  assumed  to  be  less  than  the  stationary 
demand  rate  at  each  location  so  that  occasional  outside  procurements  are 
necessary.  The  objectives  of  this  paper  are  to  develop  a  cost  model  of 
this  system  given  that  each  location  follows  a  continuous  review  procurement 
policy,  and  to  incorporate  the  single  item,  single  location  solution  method 
of  Reference  6  into  an  iterative  algorithm  which  determines  the  policy 
parameter  values  at  each  location. 


1.  INTRODUCTION 


Inventory  systems  with  returns  are  systems  in  which  there  are  units 
returned  in  a  repairable  state,  as  well  as  demands  for  units  in  a  service¬ 
able  state,  where  the  return  and  demand  processes  are  assumed  to  be 
independent.  A  discussion  of  the  contexts  in  which  this  problem  arises 
as  well  as  an  analysis  of  a  single  item,  single  location  system  are  pre¬ 
sented  in  Reference  6.  The  present  paper  extends  those  results  to  the 
management  of  a  single-item,  multi-echelon  inventory  system  consisting  of 
a  warehouse  (the  upper  echelon)  that  supports  N(N  _>  1)  retailers  (the 
lower  echelon).  In  particular  our  objective  is  to  develop  both  a  model 
for  this  problem  and  an  algorithm  for  finding  values  for  the  appro¬ 
priate  policy  variables. 

Only  one  previous  paper  has  dealt  with  the  problem  of  inventory  sys¬ 
tems  with  returns  in  a  multi -location  system.  Hoadley  and  Heyman  [3]  exam¬ 
ined  a  two-echelon  inventory  system  with  outside  procurement,  returns, 
disposals,  and  transshipment;  but  their  model  is  a  one-period  model,  and  all 
of  the  mentioned  transactions  are  assumed  to  occur  instantaneously.  Their 
approach  is  not  easily  adaptable  to  the  case  we  will  examine  in  which  repair 
and  procurement  lead  times  are  non- zero  and  the  planning  horizon  is  infinite. 

The  inventory  system  we  will  study  is  described  in  detail  in  Section  2. 
In  Section  3,  we  show  how  the  single-location  solution  method  presented  in 
Reference  6  can  be  used  to  solve  for  the  policy  parameter  values  of  a  par¬ 
ticular  cost  model  for  this  multi-echelon  system.  In  Section  4  we  conclude 
with  some  final  comments  and  suggestions  for  future  research. 
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2.  DESCRIPTION  OF  THE  SYSTEM 

As  we  have  stated  the  inventory  system  we  will  study  is  a  two-echelon 
system.  The  upper  echelon  consists  of  a  warehouse  having  both  repair  and 
storage  facilities  that  support  the  N  lower  echelon  retailers.  The  re¬ 
tailers  only  have  storage  facilities. 

All  primary  customer  demands  and  returns  are  assumed  to  occur  only  at 
the  retailers.  We  also  assume  that  all  customer  demands  not  immediately 
satisfied  are  backordered,  the  demand  and  return  processes  are  mutually 
independent  Poisson  processes,  and  lateral  resupply  is  not  allowed  between 
retailers . 

Let 


X..  =  customer  demand  rate  at  retailer  j  (j  =  1,...,N), 

Yj  =  customer  return  rate  at  retailer  j  (j  =  1,...,N), 

T^  =  constant  transportation  time  between  the  warehouse  and  a  retailer,  and 
=  constant  procurement  lead  time  for  the  warehouse  from  an  outside 
source . 

The  assumptions  that  transportation  times  are  identical  between  the  ware¬ 
house  and  any  of  the  retailers,  and  that  customer  demands  and  returns  occur 
only  at  the  retailers  are  made  for  notational  simplicity  only.  It  will 
be  apparent  that  relaxing  these  assumptions  poses  no  additional  problems . 

Recall  that  repair  facilities  exist  only  at  the  upper  echelon.  Consequent¬ 
ly,  we  assume  that  a  returned  repairable  unit  to  a  retailer  is  immediately 
sent  to  the  warehouse  from  the  retailer  and  need  not  go  back  to  that  same  re¬ 
tailer  after  it  is  repaired.  We  also  assume  that  the  repair  process  at  the 
warehouse  operates  as  a  first-come,  first-served  queueing  system. 

Since  transportation  times  are  assumed  to  be  constant,  returns  of  repair- 

N 

able  units  to  the  warehouse  occur  as  a  Poisson  process  with  rate  Yn  =  I  Yj • 

j=l  1 


Therefore,  it  is  equivalent,  and  more  convenient,  to  think  of  returns 
occurring  only  to  the  warehouse,  and  as  a  Poisson  process  with  rate  y. . 

We  assume  that  retailer  j  uses  an  (S^-l,  S_. )  continuous  review 
ordering  policy,  i.e.  a  constant  inventory  position  (net  inventory  plus 
on-order)  of  is  maintained.  This  implies  that  retailer  j  immedi¬ 

ately  orders  one  unit  from  the  warehouse  every  time  a  customer  demand 
occurs  at  that  retailer.  Since  each  order  placed  at  a  retailer  also  re¬ 
sults  in  a  demand  being  placed  upon  the  warehouse ,  demands  on  the  warehouse 

N 

occur  as  a  Poisson  process  with  rate  A 


=  Uy 
j=l  3 


[Note  the  importance  of  the  assumption  of  following  an  (S^-l,  S^) 
policy  at  retailer  j.  If  the  retailers  followed  (Q,r)  ordering  poli¬ 
cies,  then  the  time  between  the  placing  of  orders  upon  the  warehouse  would 
not  necessarily  be  exponential,  nor  would  the  orders  necessarily  be  for  in¬ 
dividual  units.  Thus  the  demand  process  at  the  warehouse  would  no  longer 
be  a  simple  Poisson  process.] 

We  assume  that  y^  >  AQ  so  that  an  occasional  outside  procurement 
is  necessary.  The  warehouse  is  assumed  to  follow  a  (Q0,r0^  P°iicY»  i.e. 
when  its  inventory  position  (net- inventory  plus  on-order  plus  in-repair) 
falls  below  rQ  +  1,  an  order  for  units  is  placed  upon  an  outside 

procurement  source.  (Note  that  a  (Q0,rQ)  policy  is  not  necessarily  opti¬ 
mal  for  reasons  discussed  in  Reference  6.  However,  we  use  it  because 
it  is  a  simple  and  commonly  used  policy.  Furthermore,  as  is  shown  in  [4], 
the  (Qq.Tq)  policy  performed  well  against  other  policies  that  were  tested.) 

Warehouse  procurement  orders  are  assumed  to  arrive  at  the  warehouse 


Tj  time  units  after  the  order  is  placed.  However,  an  order  placed  by  a 
retailer  upon  the  warehouse  does  not  necessarily  arrive  at  the  retailer 
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T  time  units  after  it  is  placed.  In  addition  to  the  transportation  time 
T  ,  there  may  be  a  delay  due  to  the  warehouse  being  out  of  serviceable 
stock.  All  demands  made  upon  the  warehouse  that  are  not  immediately  sat¬ 
isfied  are  backordered. 

A  schematic  representation  of  this  system,  is  given  by  Figure  1. 


\Ix 


Warehouse 
Return  rate 
Demand  rate 

Lead  Time  T2 

(Qo’ro>  P°licy 


Figure  1 .  A  Schematic  Representation  of  the  Multi-Echelon  System 
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Finally,  let  the  system  cost  parameters  be  as  follows: 
hQ  =  the  holding  cost  at  the  warehouse  ($/unit-year), 
hj  =  the  holding  cost  at  retailer  j  ($/unit-year)  (j  =  1,...,N), 

=  the  backorder  cost  at  retailer  j  ($/unit-year)  (j  =  1,...,N), 
and  A  =  the  fixed  warehouse  procurement  order  cost  ($/order). 

Given  values  of  h^  (j  =  ir^  (j  =  1,...,N),  and  A,  all 

assumed  to  be  non-negative,  the  problem  is  to  determine  values  for  Qq,  r^, 
and  Sj  (j  =  1,...,N)  that  will  minimize  the  expected  annual  sum  of  order¬ 
ing,  holding,  and  backorder  costs  over  all  locations. 


3.  SOLUTION  METHOD 

a)  Motivation  For  The  Algorithm 

We  begin  by  summarizing  the  results  for  the  single-item,  single- location 
problem  with  constant  procurement  lead  time  t,  Poisson  demand  rate  X, 
and  Poisson  return  rate  y  (y  <  X)  that  are  developed  in  reference  6. 


Let 


and 

x 

4>(x)  =  /  <|>(t)dt. 

t=-» 


Given  ordering,  holding,  and  backorder  costs  of  A,  h,  and  ir,  respec¬ 
tively,  and  the  fact  that  the  stationary  distribution  of  net  inventory  can 

be  closely  approximated  by  a  normal  distribution  (as  shown  in  reference  6), 
the  optimal  procurement  lot-size,  Q*,  is  the  Q  that  satisfies 
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a  =  (n+h)(j>(4> 

d  =  — ^—5-  -  •—  +  Var[lim(R(t) )]  +  (A+y)t,  and 
(A-y) 

R(t)  =  the  number  of  units  in  repair  at  time  t. 


The  optimal  reorder  point. 


satisfies 


IQ*)  ,d  JL_)  -  Iq*  -  c 


where 

c  =  -r—  +  4  +  E[lim  R(t)3  -  (X-Y)t. 

X-Y  2 

1  t-x® 

For  the  multi-echelon  problem  note  that  we  do  not  have  an  explicit 
value  for  n^,  the  warehouse  backorder  cost.  Given  the  interactions  be¬ 
tween  the  two  echelons  of  our  inventory  system,  the  cost  of  a  backorder  at 
the  warehouse  is  an  imputed  one.  It  is  measured  by  the  effect  of  a  back¬ 
order  at  the  upper  echelon  upon  the  expected  performance  at  the  lower 
echelon . 

The  optimal  stock  level  at  retailer  j,  S$  »  is  a  function  of  the 
procurement  resupply  time,  that  is, the  expected  time  from  the  placement  to 
receipt  of  an  order  by  a  retailer.  This  procurement  resupply  time  is  then 
the  transportation  time  T  ,  plus  the  expected  delay  due  to  the  warehouse 
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being  out  of  serviceable  stock.  Clearly,  costs  at  the  lower  echelon  can 
be  reduced  by  reducing  the  expected  resupply  time.  This  can  only  be  ac¬ 
complished  by  decreasing  the  expected  warehouse  backorders  at  a  random 
point  in  time,  which  is  achieved  by  increasing  Qq  and  rQ  (or  both)  at  the 
warehouse.  This,  in  turn,  raises  holding  costs  at  the  warehouse.  Thus 
a  tradeoff  exists  between  costs  at  the  upper  echelon  and  costs  at  the  lower 
echelon.  We  now  develop  an  iterative  algorithm  based  on  this  trade-off 
which  alternates  between  finding  stock  levels  for  the  upper  and  lower  eche¬ 
lons. 

Suppose  we  are  given  a  backorder  cost  for  the  warehouse,  Then 

we  can  use  (1)  and  (2)  to  find  optimal  values  for  the  parameters  Qq  and 
rQ.  These  determine  a  'performance  level'  B,  where 

B  =  the  expected  backorders  at  the  warehouse  at  a  random  point  in  time. 

2 

Letting  p  and  o  be  the  mean  and  variance,  respectively,  of  the  normal 
approximation  to  the  stationary  distribution  of  net  inventory  at  the  ware¬ 
house,  we  have 


B  =  <?<}>(— )  - 

o  o 

Then  the  expected  resupply  time  for  a  retailer  is 


(3) 


T  =  T  +  B/X  , 


It  is  shown  in  eference  4  that  the  normal  distribution  improves  as  an 
approximation  to  the  steady-state  distribution  of  net  inventory  as  ^ 
decreases  and  At  increases.  However,  it  is  a  good  approximation  in  the 
tail  of  the  distribution, which  is  all  that  we  require,  in  all  cases. 
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since  the  expected  delay  time  per  demand  is  the  expected  number  of  delays 
at  a  random  point  in  time  divided  by  the  demand  rate.  This  is  a  direct 
application  of  Little's  Formula  L  =  AW.  Then,  using  Palm's  Theorem  [1] 
as  an  approximation,  we  assume  the  number  of  units  in  resupply  at  retailer 
j  to  be  Poisson  distributed  with  mean  A_.T.  (Palm's  Theorem  requires 
the  independence  of  resupply  times ,  making  this  system  analagous  to  an 
M/G/00  queue.  Resupply  times  in  our  system  are  not  independent;  consider, 
for  example,  a  demand  by  a  retailer  which  cannot  be  immediately  filled  by 
the  warehouse.  Then  it  is  more  likely  that  the  next  demand  placed  by  a 
retailer  upon  the  warehouse  also  experiences  a  delay  than  if  the  preceding 
order  had  been  immediately  satisfied.) 

The  use  of  Palm's  Theorem  to  approximate  the  distribution  of  the 
number  of  units  in  resupply  at  retailer  j  with  a  Poisson  distribution 
is  discussed  later.  (As  we  will  see,  this  approximation  is  good  when  the 
expected  warehouse  net  inventory  is  non-negative . ) 

Once  we  know  the  value  of  T  and  have  the  form  (approximately)  of  the 


distribution  of  the  number  of  units  on-order  by  retailer  j  we  can  solve 


N  independent  subproblems  to  obtain  the  optimal  value  for  S ^ .  The  sub- 
problem  at  retailer  j  consists  of  finding  the  optimal  stock  level  S?, 
assuming  a  constant  procurement  resupply  time  of  T,  where  T  is  given 


by  (3).  This  is  accomplished  using  Lemma  1. 


Lemma  1.  Suppose  the  procurement  lead  time  is  a  constant  T  and  demand 
is  Poisson  distributed  with  rate  A _. .  Then  the  optimal  value  S*  for  an 
(Sj-1,  Sj)  policy  is  the  largest  integer  S^  such  that 
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(4) 


h. 

P(S . ;  X.T)  >  v  1 
—  -|  -i  tt  .  +h . 

J  J  3  3 


where 


P(x,u)  =  l  p(r;p) 
r=x 


and 


(  \  -  "4  4 

p(r;u)  -  e  — p 


The  proof  of  Lemma  1  can  be  found  on  page  204  of  Reference  2. 

Let  K  (S  ,T)  be  the  expected  annual  holding  and  backorder  costs  at 

j  j 

retailer  j  when  the  constant  inventory  position  is  S_.  and  the  constant 
procurement  lead  time  is  T.  As  can  be  shown  (see  Reference  2) 


(5)  K  •  ( S  . ,  T )  =  (ir.+h.)[X.TP(S.-l;  X.T)  -  S.P(S.;  X.T)]  +  h.[S  -  X  Tj. 

^  3  3^3  3  3  333  333 


For  a  fixed  value  of  T  (and  therefore  of  B)  we  define  the  minimum  total 

£ 

expected  costs  at  the  lower  echelon,  K  (B),  as 


(6) 


N 

K  (B)  =  l  K.(S*,T), 

j=l  :  1 


where  K.(*,*)  is  given  by  (5),  T  is  given  by  (3),  and  S!$  satisfies 


(4).  Given  a  current  value  of  B  =  b. 


dK*(B) 


dB 


3 


will  be  our  current 


B=b 


estimate  of  ,  since  it  measures  the  marginal  effect  of  a  warehouse  back¬ 
order  on  the  expected  total  lower  echelon  cost.  It  is  easy  to  show  that 


(7) 


dK*(B) 


-4^-  =  r-  1  [(*.+  h.)X.P(S*;  X.T)  -  h.X  ]. 

dB  *0  -5=1  J  3  3“  3  3  3  3 


b)  The  Algorithm 

Before  presenting  its  details  we  will  give  a  brief  overview  of  the 


I 
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algorithm's  steps.  The  algorithm  begins  assuming  a  value  for  ttq  is 

available.  Given  this  value,  the  optimal  values  for  Qq  and  rQ  for  this 

single  location  problem  with  returns  are  found  using  (1)  and  (2).  This 

solution  determines  a  performance  level  B  =  b,  and  therefore  an  expected 

resupply  time  of  T  for  the  lower  echelon.  Given  this  value  of  T,  op- 

£ 

timal  stock  levels  are  determined  for  each  retailer,  and  we  can  find  K  (b), 
the  total  expected  cost  for  the  lower  echelon.  A  new  estimate  is  obtained 

),  and  the  cycle  repeats  until  convergence  occurs. 


for  *0  (tt0 


dK*(B) 


dB 


B=b 


Next,  let  K  (B)  be  the  minimum  total  of  expected  warehouse  order¬ 
ing  and  holding  costs  that  can  be  achieved  given  that  the  expected  ware¬ 
house  backorders  is  B. 

We  state  two  lemmas  without  proof. 


Lemma  2.  K  (B)  is  convex  decreasing  in  B. 


Lemma  3.  Let  T  be  a  constant  resupply  time.  If  the  optimal  stock  levels 


S^.  (j  =  1,...,N)  are  continuous  rather  than  integer- valued,  then  K  (B) 
is  a  concave  increasing  function  of  B,  where  B  =  XQ(T-T^). 


The  proof  of  these  lemmas,  which  involve  nothing  more  than  applying  the  chain 
rule  to  take  derivatives,  may  be  found  in  Reference  4. 

t  £  U 

Figure  2  represents  a  typical  graphing  of  K  (B)  and  K  (B)  against 

B.  The  objective  is  to  minimize  K  (B)  +  K  (B)  with  respect  to  B.  We 

observed  in  all  test  cases  that  under  the  conditions  of  Lemma  3,  KU(B)  + 

a 

K  (B)  was  convex  in  B.  Therefore,  the  minimum  cost  will  occur  where 

dK*(B)  _  dKU(B) 

dB  "  dB 
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(4) 


P(S . ;  X.T)  > 

-  1  D  ^j+hj 


where 


and 


P(x,p)  =  l  p(r;w) 


r=x 

-  "4  4 


p(r;p)  =  e 


The  proof  of  Lemma  1  can  be  found  on  page  204  of  Reference  2. 

Let  K  (S  ,T)  be  the  expected  annual  holding  and  backorder  costs  at 

j  j 

retailer  j  when  the  constant  inventory  position  is  and  the  constant 

procurement  lead  time  is  T.  As  can  be  shown  (see  Reference  2) 


(5)  Kj(Sj  ,T)  =  (ir  +h  )[A^TP(S^-1;  X..T)  -  S^P^CS^ ;  X..T)]  +  1k[S^-  X_.T]. 


For  a  fixed  value  of  T  (and  therefore  of  B)  we  define  the  minimum  total 

jj 

expected  costs  at  the  lower  echelon,  K  (B),  as 


(6) 


N 

K  (B)  =  l  K.(S*,T), 
j=l  3  3 


where  Kj(*,’)  is  given  by  (5),  T  is  given  by  (3),  and  S$  satisfies 
(4).  Given  a  current  value  of  B  =  b. 


dK*(B) 


dB 


will  be  our  current 


B=b 


estimate  of  since  it  measures  the  marginal  effect  of  a  warehouse  back¬ 

order  on  the  expected  total  lower  echelon  cost.  It  is  easy  to  show  that 


(7) 


dK*(B) 

dB 


l  N 

r~  l  [(*.+  h.)X.P(S*;  X.T)  -  h  .X .]. 


o  j=i 


D  1 


:  i 


b)  The  Algorithm 


Before  presenting  its  details  we  will  give  a  brief  overview  of  the 


Figure  2.  Minimum  Upper  and  Lower  Echelon  Cost  Functions  Vs.  B 


As  was  shown,  it  is  easy  to  find  the  optimal  stock  levels  S*,  using 


Lemma  1.  Using  these  stock  levels  we  find  K  (B)  for  a  fixed  value  of 


B.  Finding  K  (B)  directly  as  a  function  of  B  is  extremely  difficult 


Thus 


B  at  which  the  two 


curves  in  Figure  2  have  tangents  with  slopes  of  opposite  sign.  The  slope 


of  K  (B)  at  this  optimal  value  of  B  is  the  imputed  warehouse  backorder 


Suppose  we  originally  had  this  value  of  n  and  solved  the  ware 


house  problem  (using  (1)  and  (2)).  The  optimal  value  of  B  for  this  problem 


would  be  the  same  value  of  B  that  minimizes  K  (B)  +  K  (B) 


The  first  few  steps  of  the  algorithm  are  illustrated  in  Figure  3.  The 


algorithm  begins  by  setting  n  =  max  (tt.).  This  is  an  upper  bound 

j=l,...,N  3 

on  the  optimal  value  of  ir  ,  since  this  value  implies  that  a  backorder 


at  the  warehouse  always  results  in  a  backorder  at  the  retailer  with  the 


largest  backorder  cost.  Then  Qn  and  r  are  found  using  this  upper  bound 


This  determines  a  value  of  B 


which  is  a  lower  bound  on  the  optimal  value  of  B  (and  therefore  of  T) 


These  computations  yield  point 


on  the  upper  echelon  cost  curve  in 
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Using  this  lower  bound  on  the  optimal  value  of  T,  we  find  a  lower 

bound  estimate  of  S*  (j  =  1,...,N)  which  determines  a  value  K  (b  ), 

_  3  o  .  ,  A 


point  ^2^  in 


Figure  3.  Next  we  set  ir. 


_  dK*(B) 


dB 


Since  K  (B) 


B=b, 


is  concave  in  B,  and  since  we  have  a  lower  bound  estimate  of  the  op¬ 
timal  B,  the  new  estimate  of  is  an  upper  bound  on  the  optimal  value 

of  but,  it  is  smaller  than  the  previous  estimate.  Using  this  new 

estimate  of  tt^,  B  will  increase  to  a  value,  say  b2,  as  a  result  of 

resolving  for  Q  and  r  using  (1)  and  (2).  These  calculations  produce 
u  u 

point  (3)  in  Figure  3.  The  procedure  continues  by  letting  T  =  T  +  -p- 

w  £  r\  0 

and  finding  K  (b2>,  which  leads  to  point  (4)  .  The  algorithm  continues 
in  this  manner  until  convergence  occurs. 


Figure  3.  First  Steps  of  Algorithm 


The  algorithm  we  have  described  can  formally  be  stated  as  follows: 


Algorithm 


max 


r  using  Equations  (1)  and  (2)  and 


Given  n, 


determine  the  corresponding  value  of  B,  say  b 


find  the  S*  using  Equation  (4) 


Step  3)  Using  these  S$,  find 


evaluated  at  B  =  b,  using 


Equation  (7) 


let  ir  assume  this  value,  and  return  to  Step 


That  the  algorithm  converges  can  be  observed  by  considering  the 


problem  as  the  following  constrained  optimization  problem 


If  we  let  n  be  the  Lagrange  multiplier  for  the  constraint  in 


problem  PI  (i.e.  it  is  the  imputed  cost  of  a  warehouse  backorder),  the 
above  problem  may  be  reformulated  as  the  following  optimization  problem 
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The  optimality  conditions  for  this  problem  are 


(8) 


d{  l  K.(S*,T)} 

--3*i  3  3 

dT 


-  = 

0  0  ’ 


(9) 


dB  0 


and 


(10) 


B  -  ^0(T-T1)  =  0. 


Since  K  (B)  wa*;  defined  to  be  equal  to  [  K.(S*,T), 


j=i 


0  d{  l  K  (S*,T)}  d{  i  K  (S*,T)} 

(11)  1  ’  .  a ,  ]»i  ^  j 


dB 


dT 


dB 


dT 


»0  •  V 


Thus  from  (9)  and  (11)  we  see  that 


dK*(B)  _  dKU(B) 


dB 


dB 


=  w 


In  Step  2)  of  the  algorithm,  condition  (10)  is  maintained,  while 
condition  (11)  is  not.  The  value  of  B,  determined  by  the  values  of 
Qq  and  rQ  in  the  previous  step,  establishes  the  value  of  T^.  However, 

£  II 

the  slopes  of  the  two  curves  K  (B)  and  K  (B)  are  not  necessarily  of 
opposite  sign  at  this  value  of  B. 


value  of 


In  Step ^3)  and  then  Step  1)  the  reverse  is  true.  In  using  the 

~ as  the  warehouse  backorder  cost,  condition  (11)  is 
B=b 


dB 
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maintained,  while  condition  (10)  is  not. 

However,  since  the  value  of  B  is  monotonically  increasing  (and 

is  monotonically  decreasing)  with  every  iteration  of  the  algorithm, 

the  value  of  B  converges  on  that  value  such  that  both  (10 )  and  (11 ) 

are  satisfied.  The  value  of  converges  to  the  optimal  value  of  the 

Lagrange  multiplier  in  Problem  (PI).  An  example  showing  the  results 

of  applying  the  algorithm  are  given  in  Table  1. 

Two  complications,  however,  occurred  in  a  few  of  the  50  test  cases 

£ 

that  were  run.  Recall,  first,  that  Lemma  3  states  that  K  (B)  is  con¬ 
cave  increasing  if  the  stock  levels  S$  are  allowed  to  vary  continuously. 
Actally,  Lemma  1  is  used  to  find  an  integer  value  of  J*.  For  fixed 

values  of  Sj  (j  *  1,...,K),  K  (B)  is  actually  locally  convex  in  B 

almost  everywhere  and  not  concave,  even  if  the  'general  shape’  of  the 

£ 

curve  is  still  concave.  A  true  representation  of  the  curve  K  (B)  is 

given  by  Figure  4.  Thus  when  restricting  ourselves  to  integer  values 

for  the  S*,  an  overestimate  rather  than  an  underestimate  of  in 

Step  3)  of  the  algorithm  is  possible.  Even  so,  convergence  to  the  correct 

value  of  7rQ  (and  of  B)  is  still  to  be  expected.  This  overestimation 

of  ^  did  occur  in  a  few  of  the  test  cases. 

In  general,  the  optimal  values  of  Qq,  rQ  and  S_.  ( j  =  1,...,N) 

were  found  after  only  three  iterations  of  the  algorithm.  This  occurred 

in  48  of  the  50  test  cases.  The  curve  K  (B)  is  very  flat  compared  to 

KU(B),  so  that  convergence  to  the  correct  value  of  nQ  occurs  quickly. 

£  11 

Recall  that  we  noted  earlier  that  K  (B)  +  K  (B)  was  convex  for  all  of 

0 

the  50  test  problems.  The  reason  this  occurred  was  because  K  (B),  al¬ 
though  concave,  is  almost  linear. 
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Table  1 


An  Example  of  Algorithm  Intermediate  Results 
N  =  5,  =  20,  Y0  =  8,  c  =  -15,  d  =  40,  A  =  20,  ^  =  1 

hj  =  1,  n..  =  10,  (j=l,...,5),  hQ  =  .6 

A1  =  2,  X2  =  3,  X3  =  4,  X4  =  5,  X5  =  6 


Iteration  it 

*o 

Q0 

ro 

B 

T 

s*(j 

=  1, 

•  *  * 

,5) 

1 

10.00 

14 

20 

.124 

1.01 

4 

5 

7 

8 

9 

2 

0.59 

21 

4 

3.708 

1.19 

5 

6 

8 

9 

11 

3 

0.45 

22 

4 

3.791 

1.19 

5 

6 

8 

9 

11 

17 


c)  Testing  the  Poisson  Assumption 

Recall  that  in  the  development  of  the  algorithm  we  approximated 

the  distribution  of  the  number  of  units  in  resupply  at  retailer 

j  (j  =  1,...,N)  with  a  Poisson  random  variable  with  mean  X^T.  To 

test  this  approximation,  we  derived  an  expression  for  R^(t),  the 

number  of  units  in  resupply  at  retailer  j,  for  a  special  case,  namely 

one  in  which  the  repair  facility  at  the  warehouse  behaves  as  an  M/D /“ 

queueing  system.  The  derivation  was  similar  to  ones  given  in  References 

5  and  7;  the  details  may  be  found  in  Reference  4. 

Using  this  exact  distribution  of  R^(t)  the  assumption  that 

R.(t)  is  Poisson  with  mean  X  t  was  tested.  As  was  found  in  Reference  7, 
J  D 

the  results  indicate  that  the  Poisson  assumption  (i.e.  the  assumption 

that  resupply  times  are  virtually  independent)  improves  as  the  expected 

warehouse  backorders,  or  the  probability  of  delay  at  the  warehouse, 

decreases.  In  particular,  the  Poisson  assumption  was  found  to  be  good 

as  long  as  the  expected  value  of  net  inventory  at  the  warehouse  at  a 

random  point  in  time  is  greater  than  zero.  This  will,  of  course,  be  the 

case  for  a  reasonably  large  ratio  of  backorder  to  holding  costs. 

A  sample  system  showing  how  the  Poisson  approximation  to  R^(t) 

improves  as  the  warehouse  'safety  stock'  increases,  is  given  in  Figures 
N 

5-7.  Let  X'  =  J  X.,  and  let  R  be  the  constant  warehouse  repair 
0  i=i  1 

i/i 

time.  In  our  example,  suppose  that  X^  =  3,  X^  =  2,  =  1,  T^  =  1, 

T2=2,  R=.6,  S^=3,  and  Qq  =  5.  The  values  of  rQ  considered  are 
4,6,  and  8.  Table  2  shows  how  some  of  the  performance  measures  change 


with  these  values  of  rQ. 


18 


Table  2 

Performance  Measures  as  rQ  Changes 


ro 

B 

T 

R^(t)  vs.  Poisson  Approximation 

4 

1.462 

1.292 

Figure  5 

6 

.732 

1.146 

Figure  6 

8 

.318 

1.064 

Figure  7 

If  the  expected  value  of  net  inventory  at  the  warehouse  is  to  be 
greater  than  zero,  when  Qq  is  fixed  at  5,  then  rQ  must  be  at  least 

7 

5  — .  When  rQ  =  6,  the  Poisson  approximation  to  the  number  of  units  in 
resupply  at  base  j,  R..(t),  is  already  quite  good.  When  rQ  =  8,  the 
approximation  is  almost  perfect. 

Another  complication  that  occurred  in  our  50  test  cases  was  due  to 
the  fact  that  the  Poisson  approximation  to  R^(t)  underestimates  the  variance 
when  B,  the  expected  number  of  warehouse  backorders,  is  large.  Thus  the 
use  of  the  Poisson  approximation  in  defining  K  |b)  by  (6)  serious¬ 
ly  underestimates  the  true  expected  total  cost  at  the  lower  echelon  when 
B  is  large.  In  a  few  of  the  test  cases,  low  values  of  rQ  and  Qq 
were  assigned  to  the  warehouse,  which  made  its  expected  inventory  posi¬ 
tion  less  than  its  expected  net  lead  time  demand .  This  occurred  because 
the  true  lower  echelon  expected  costs  for  relatively  large  values  of 
B  were  significantly  underestimated. 

One  can  make  a  simple  alteration  in  the  algorithm  to  alleviate  this 
difficulty.  It  is  reasonable  to  require  that  the  expected  value  of  net 


probability 


.20 


.15 


.10 


.0 


5 


inventory  at  the  warehouse  at  a  random  point  in  time  be  greater  than 
zero.  From  Equation  (16)  of  Reference  6,  we  see  that  h^  assures 

us  that  this  condition  will  hold.  Thus,  we  can  alter  Step  3)  of  the 
algorithm  so  that  tt  is  set  equal  to  the  maximum  of  the  expression 
given  in  Step  0)  of  the  algorithm  and  h^.  This  assures  us  of  having 
reasonable  warehouse  stock  levels  and  of  having  the  Poisson  distribution 
be  a  good  approximation  to  the  distribution  of  the  number  of  units  in 
resupply  at  a  retailer. 

4.  CONCLUDING  COMMENTS 

We  have  shown  how  the  single-location  solution  method  of  [6]  may 
be  incorporated  into  an  iterative  algorithm  for  setting  stock  levels  in 
the  single-item,  multi-echelon  inventory  problem  with  returns.  In  [4] 
it  is  shown  that  simple  alterations  allow  this  algorithm  to  set 
stock  levels  in  a  multi-echelon  system  in  which  repair  of  returned  units 
may  also  occur  at  the  retailers.  A  necessary  condition  is  that  there 
be  a  large  number  of  retailers,  so  that  the  return  process  at  the  ware¬ 
house  is  well  approximated  by  a  Poisson  process. 

The  algorithm  presented  may  also  be  imbedded  in  a  larger  algorithm 
which  sets  stock  levels  in  an  M-echelon  inventory  system  with  returns. 

The  only  requirement  would  be  the  use  of  an  (S-l,  S)  procurement  policy 
at  each  of  the  lower  M-l  echelon  locations. 

We  conclude  with  a  discussion  of  possible  directions  for  future 
research  on  the  problem  of  inventory  systems  with  returns.  Though  the 
inventory  models  considered  in  this  paper  and  in  Reference  6  contain  less 
restrictive  assumptions  than  those  of  the  previous  literature,  many  im¬ 
portant  extensions  still  need  to  be  examined.  For  the  single- location 
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problem,  these  include  non- stationary  return  and  demand  processes,  a 
finite-horizon  periodic  review  model,  and  a  multi-item  model  in  which 
repairable  units  of  different  types  'compete'  for  the  limited  repair 
capacity.  For  the  multi-echelon  case,  an  analysis  should  be  conducted 
of  the  warehouse  demand  process  when  the  retailers  follow  different 
(Q,r)  procurement  policies,  and  an  algorithm  needs  to  be  developed 
for  determining  procurement  policies  for  the  warehouse  and  the  retailers 
when  all  locations  use  a  periodic  review  policy. 


24 


BIBLIOGRAPHY 

[1]  G.J.  Feeney  and  C.C.  Sherbrooke,  "The  (S-l,  S)  Inventory  Policy 

Under  Compound  Poisson  Demand,"  Management  Science  12,  p.  391- 
411  (1966). 

[2]  G.  Hadley  and  T.M.  Whitin,  Analysis  of  Inventory  Systems,  Prentice- 

Hall,  New  Jersey  (1963). 

[3]  B.  Hoadley  and  D.P.  Heyman,  "A  Two-Echelon  Inventory  Model  with 

Purchases,  Disposition,  Shipments,  Returns,  and  Transhipments," 
Naval  Research  Logistics  Quarterly  24,  p.  1-19  (1977). 

[4]  M.H.  Isaac,  "An  Analysis  of  Inventory  Systems  With  Returns," 

Unpublished  Ph.D.  dissertation.  School  of  Operations  Research  and 
Industrial  Engineering,  Cornell  University  (1979). 

[5]  J.A.  Muckstadt,  "Analysis  of  a  Two-Echelon  Inventory  System  In 

Which  All  Locations  Follow  Continuous  Review  (S,s)  Policies," 
Technical  Report  No.  337,  School  of  Operations  Research  and 
Industrial  Engineering,  Cornell  University,  Ithaca,  N.Y.  (1977). 

[6]  J.A.  Muckstadt  and  M.H.  Isaac,  "An  Analysis  of  a  Single  Item  In¬ 

ventory  System  With  Returns:  The  Single  Echelon  Case,"  submitted 
to  Naval  Research  Logistics  Quarterly. 

[7]  K.  Shanker,  "An  Analysis  of  a  Two-Echelon  Inventory  System  For 

Recoverable  Items,"  Technical  Report  No.  341,  School  of  Operations 
Research  and  Industrial  Engineering,  Cornell  University,  Ithaca, 
N.Y.  (1977). 


UNCLASSIFIED 


security  classification  of  this  face  rwsm  o«<  e»i«o 


I  REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

T-Sl-oWf  NUMECN - 

432  ' 

t.  GOVT  ACCESSION  NO. 

1.  RECIPIENT’S  CATALOG  NUMECR 

4.  TITLE  (end  Sukttllo) 

AN  ANALYSIS  OF  SINGLE  ITEM  INVENTORY  SYSTEM 

WITH  RETURNS:  THE  MULT I -ECHELON  CASE 

TECHNICAL  REPORT 

S.  PERFORMING  ORO.  REPORT  NUMBER 

7.  AUTHOR/ «J 

JOHN  A.  MUCKSTADT 

MICHAEL  H.  ISAAC 

S.  CONTRACT  OR  GRANT  NUMBER/*/ 

N00014-75-C-1172  J 

A  PERFORMING  ORGANIZATION  NAME  ANO  AOONESS  J 

SCHOOL  OF  OPERATIONS  RESEARCH  AND  INDUSTRIAL 
ENGINEERING 

CORNELL  UNIVERSITY,  ITHACA,  NY  14853 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  B  PORK  UNIT  NUMBERS 

NR  042-335 

It.  CONTROLLING  OFFICE  NAME  ANO  AOORCSS 

STATISTICS  AND  PROBABILITY  PROGRAM 

OFFICE  OF  NAVAL  RESEARCH  (Code  436) 

ARLINGTON,  VA  20360  ~ 

IS.  REPORT  DATE 

JULY  1979 

IS.  NUMBER  OF  PAOES 

24 

16.  MONITORING  AGENCY  name  k  ADDRESS/!/  dllloiunt  Item  Cenuelllnc  Olll cm) 

IS.  SECURITY  CLASS,  ft  thf  rmpmtt) 

UNCLASSIFIED 

l4«L  DECLASSIFICATION/  DOWNGRADING 
SCHEDULE 

IS.  OISTNISUTION  STATEMENT  (ol  Mila  XaFWIJ 


APPROVED  FOR  PUBLIC  RELEASE:  DISTRIBUTION  UNLIMITED. 

17.  OISTNISUTION  STATEMENT  (ol  the  okotroet  entered  In  Block  10.  II  dlllotonl  trom  Retort) 


is.  supplementary  notes 

PRESENTED  AT  ORSA/TIMS  CONFERENCE,  NEW  ORLEANS,  MAY  1979 

•t.  NET  WOROS  (Continue  on  toeoroo  oldo  II  no cooomr  end  idontltr  kr  klock  mmkor) 

INVENTORY  THEORY,  REPAIRABLE  ITEM,  OPTIMAL  POLICY,  MULTI-ECHELON 

ASST  NACT  (Continue  on  tooo too  cldc  II  nocooomr  end  I  dent  I  If  kf  Clock  mmtkor) 

We  consider  the  inventory  control  of  a  single  item  in  a  two  echelon  system 
in  which  a  warehouse  (the  upper  echelon)  supports  N(N  Tpl )  retailers 
(the  lower  echelon).  Customers  return  units  in  a  repairable  state  as  well 
as  demand  units  in  a  serviceable  state  at  the  retailer  level  only.  The 
return  and  demand  processes  are  assumed  to  be  mutually  independent .  The  ! 
stationary  return  rate  is  assumed  to  be  less  than  the  stationary  demand 
rate  at  each  location  so  that  occasional  outside  procurements  are  necessary, 


00  ,ET„  1473 


COITION  OF  I  NOV  <*  IS  OMOLCTE 
S/N  0102-  LF-014.  6401 


_ _ _ ^ 

SECURITY  CLASSIFICATION  OF  THIS  e  cat  ,'Am  O.l*  Bntorod) 


UNCLASSIFIED 


4 


UNCLASSIFIED 


•CCUMTV  CLASSIFICATION  or  THIS  r»a«  r*S*>  Omm  bims 


Abstract  -  continued 

"^The  objectives  of  this  paper  are  to  develop  a  cost  model  of  this  system 
given  that  each  location  follows  a  continuous  review  procurement  policy, 
and  to  incorporate  the  single  item,  single  location  solution  method 
^-<jf  Reference.  G  into  an  iterative  algorithm  which  determines  the  policy 
parameter  values  at  each 'location. 


\ 


>  'i  o;o;.  ju* *«o: 


UNCLASSIFTrn 


ICCjNlT''  Cl»Mi»'CaTiOn  sr  '-.I  3m •  Cni*r*«l 


